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ABSTRACT 

A preliminary study was made to define problem areas and potential advantages of a 
radiator for extraterrestrial applications in which a heat exchange liquid is circulated 
through a capillary medium. An evaluation of surface energy, heat transfer, hydraulic, 
geometric, and meteoroid penetration factors in a capillary radiator operating under 
extraterrestrial conditions is presented to indicate feasibility in a potential application. 
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A RADIATOR CONCEPT BASED ON CAPILLARY RETENTION OF LIQUID 

by Alex Vary 

Lewis Research Center 


SUMMARY 

A radiator concept in which a heat exchange liquid is circulated in a capillary me- 
dium contained in heat radiating panels is described. It is proposed that this concept 
may be applied to space systems where the possibility of meteoroid punctures exists. 

An argument based on well-known surface energy phenomena is presented to show that 
capillary retention of liquid at puncture sites can be expected to halt liquid leakage under 
certain conditions. Liquid escape by the alternative mode of vaporization can be sub- 
stantially reduced by employing liquids having low vapor pressures and by operating the 
radiator at sufficiently low temperatures. 

The feasibility of a radiator based on the capillary concept operating under extrater- 
restrial conditions is examined in terms of surface energy, heat transfer, fluid dynamic, 
geometric, and meteoroid penetration factors. Estimates of the significant character- 
istics of a hypothetical capillary radiator indicate that the concept is potentially com- 
petitive with other concepts being considered for space application. However, a firm 
conclusion concerning feasibility must be deferred until there is experimental evidence 
of the viability of liquid-filled capillary media subjected to the conditions of an actual 
space application. Problems associated with the concept are defined, and resolutive 
investigations are suggested. 

Because of the preliminary nature of this study, predictions of the potential advan- 
tages of using a capillary radiator are speculative. Nevertheless, some possible advan- 
tages include lower specific weight, lower temperature, higher power cycle efficiency, 
and less restriction on the radiator size relative to the conventional fin -tube concept. 


INTRODUCTION 

Liquid- or vapor -filled radiators that are being considered for space applications 
should incorporate features that provide adequate protection against the deleterious 



effects of meteoroid impacts. Present design concepts include the use of armor, 
"bumper” plates, or reinforced tube walls to reduce vulnerability to meteoroid punc- 
tures. In large space powerplants the radiator weight with respect to useful power may 
be significantly higher than that of any other major component in interplanetary vehicles 
(ref. 1). Recent efforts in radiator design therefore have been directed toward reducing 
not only the vulnerability to meteoroid punctures but also the weight penalty arising from 
meteoroid protection requirements (refs. 2 and 3). 

The radiator concept that is currently contemplated for space powerplants is based 
on the circulation of heat exchange fluid in a fin-and-tube structure protected by rein- 
forcement or armor (refs. 4 and 5). Other radiator concepts that have been proposed are 
described in references 2, 3, and 6. The concept proposed herein is based on a capillary 
structure that is inherently vulnerable to meteoroid puncture but relies on capillary re- 
tention of the heat exchange liquid at anticipated puncture sites. 

In the proposed "capillary radiator" concept, the prime heat radiating surfaces con- 
sist of thin panels containing a capillary medium within which the heat exchange liquid 
circulates. When the panels are damaged by meteoroids, surface energy phenomena 
should prevent leakage if certain conditions are satisfied. Liquid loss which would then 
occur by evaporation can be reduced by selecting liquids having low vapor pressures and 
by operating the radiator at temperatures that produce low vaporization rates. 

The purpose of this report is to indicate the feasibility of a liquid -filled, capillary 
radiator that will experience numerous punctures by meteoroids. Possible advantages 
and problems associated with a capillary radiator are also discussed. 


FUNDAMENTAL CONSIDERATIONS 
General Description 

The capillary radiator consists of a number of thin extensive panels containing a 
capillary medium saturated with a heat exchange liquid. The panels are attached to 
pipes that convey the liquid between the radiator and a heat exchanger by means of a 
suitable pumping and distribution scheme. 

Internally the panels may consist of layers of sheet metal separated by spacers or 
of foam metal sandwiched between surface layers of sheet metal, as illustrated in fig- 
ure 1. There are many other equivalent modes of capillary construction; for example, 
laminations of corrugated sheet metal, laminations of wire mesh screens, metal fiber 
wicks, metallic sponges, and so forth. The crucial feature of a suitable capillary me- 
dium would be an intercommunicating cellular substructure. This "connected" cell 
feature is necessary to provide flow channels around discontinuities and obstructions 
that arise from meteoroid impacts. 
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Figure 1. - Segment of typical capillary radiator. 


Analysis of Capillary Retention 

The purpose of this section is to identify surface energy factors that influence capil- 
lary retention of liquid flowing past punctures in a capillary medium. It will be shown 
that the design of a capillary radiator depends on the specification of a capillary dimen- 
sion and the imposition of certain conditions on the capillary medium and circulatory 
system. 

Basic phenomena and prerequisite con ditions. - To obtain a general expression for 
specifying the capillary dimension, it is necessary to consider the surface tension effects 
that are discussed hereinafter. 

Throughout the following discussion, it is assumed that the equilibrium contact 
angle lies between 0° and 180°. A crucial factor to the capillary radiator concept 
is the prerequisite condition that the cohesive forces in the liquid be greater than the 
adhesive forces with respect to the container material (i. e. , no spontaneous spreading). 
The following discussion pertains to a zero- or low -gravity vacuum environment. 

Figure 2 shows a cross section of the liquid near the end of a cylindrical capillary 
tube with diameter 6 g . The illustration represents a M clean break” or right truncation 
of one element of a radiator consisting of an array of tubular channels. Under zero 
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(a) Meniscus inside capillary. (b) Meniscus at edge. 

Figure 2. - Cross section of liquid-filled capillary tube in zero gravity. 


gravity, the curvature of the meniscus depends solely on surface tension forces. 

If the liquid is nonwetting, a convex meniscus as depicted in figure 2(a) will form 
upon exposure of the liquid surface and the contact angle will be greater than 90°. If 
the liquid surface is displaced (as, for example, by increasing the liquid inventory) 
so that the meniscus is at the lip of the tube (fig. 2(b)), then surface tension effects 
restrain the liquid from climbing over the edge of the tube. This restraint exists as 
long as the angle at the edge 9 does not exceed the contact angle increased by 90°. 
A familiar example of the ’’edge effect" (see ref. 7) is the surface tension phenomenon 
that allows overfilling of a tumbler of water. A discontinuity such as an edge will bar 
movement of a liquid surface and allow the meniscus curvature to change with the liquid 
pressure. 

The angle 6 is measured between the surface and tangent to the meniscus at the 
edge as indicated in figure 2(b); 6 is equal to or greater than 9 Q the equilibrium con- 
tact angle. Although the contact angle is an invariant depending on degree of wetting 
(ref. 8), the edge angle 6 in figure 2(b) can be varied by varying the liquid inventory 
or pressure. The maximum value of 9 is 9 Q + 90° as noted in the previous paragraph. 

A "contact angle hysteresis" (refs. 7 and 8) occurs with wetting as well as non- 
wetting liquids whenever an interface tends to move. Because of the edge and hysteresis 
effects, the liquid pressure may be different from the equilibrium pressure established 
by the usual contact angle. The surfaces and edges needed to exploit these effects are 
assumed to exist in the capillary matrix. 

The previously described phenomena may be applied to liquid retention at meteoroid 
puncture sites in a capillary medium. For example, when a puncture occurs and a quan- 
tity of liquid is lost, if this quantity is not replaced, the liquid surface will reside more 
or less within the capillary. Then, as in the case illustrated in figure 2(b), a stable 
meniscus will form for a range of pressures provided that suitable surfaces and edges 
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are available. The edge and hysteresis effects will tend to bar further liquid egress 
after the initial liquid loss in the absence of large pressure changes. 

The pressure that can be supported within the liquid by surface tension forces de- 
pends upon the capillary dimension. For a tube with the meniscus at the edge as in 
figure 2(b), the capillary diameter for which the force of surface tension balances the 
force of liquid pressure is given by the following equation: 

P = — (-cos 9) (1) 

6 e 

where S g is the equilibrium diameter associated with the surface tension a, pressure P 
of the liquid, and the variable angle 9. (Symbols are defined in appendix A. ) Equa- 
tion (1) applies only to cases where 9 is greater than or equal to the equilibrium con- 
tact angle (i. e. , 0 c < 0 < 6> c + 90°). The external pressure is assumed zero throughout 
this discussion. 

Given 6 , a, and -cos 9, equation (1) sets the maximum pressure at any puncture. 
Note that P may be zero or negative as when the break meniscus curvature is concave. 
Note also that for circular tubes as in figure 2(b) the maximum positive pressure occurs 
when the meniscus is a convex hemisphere. For either wetting or nonwetting liquids, 

9 will be between 0° and 180°. 
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(a) Puncture with d > 6^. (b) Puncture with d < 6g. 

Figure 3. - Hypothetical circular puncture in capillary radiator consisting of two parallel plates. 


Selection of the capillary structure. - It is instructive to consider menisci which 
may arise at hypothetical puncture sites not having circular tubular openings. Consider, 
for example, an elemental radiator formed of two closely spaced parallel plates with 
spacing 5 T e and containing liquid at pressure P. In figure 3 a cross section of a right- 
circular hole with diameter d is shown under conditions where the interface is at the 
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break edge. If d is greater than 5^, the meniscus will be as shown in figure 3(a). 

For equilibrium 6' is determined in accordance with the following equation: 

6' = — (-cos 6 ) (2) 

e P 

In general, for circular holes having d > 6’ as in figure 3(a) the critical dimension is 5’ , 

G 6 

the spacing between plates. Figure 3(b) represents the alternative case where d<6 e 
and the "cylindrical” meniscus of figure 3(a) merges to form two "circular" menisci. 

For holes having d < 5^ as in figure 3(b), the critical dimension is d, the diameter of 
the hole. Should the hole diameters d invariably be less than 5 g as defined by equa- 
tion (1) which is applicable to figure 3(b), then the radiator could consist of parallel 
plates having any convenient spacing. Because punctures with d > 6 g will invariably 
occur, 6^ is bound to be the critical dimension in general. 

Unlike the previous idealized case the perimeters of actual punctures will usually 
be noncircular and nonplanar. For example, consider an elemental radiator formed of 
two parallel plates having an irregularly shaped hole in the form of a right parallelepiped 
as illustrated in figure 4. The hole "diameter" exceeds the plate spacing 5’ in this 
case. 



Figure 4. - Irregular puncture in capillary radiator consisting of two 
parallel plates. 
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Because of impact distortion, the plate spacing near holes will not be uniform, and 
generally the meniscus perimeter will be nonplanar. At an edge the force balance for a 
typical differential element of the meniscus (see fig. 4) with dimensions d l and 5 TT 
will vary with 5 f , the local plate spacing; so that, 

PS 1 d l = 2a(-cos 0) d l + 2a(-cos 0 T )5 ,f 

Because the angle 6 1 relative to contiguous elements of the meniscus is 90°, as indi- 
cated in figure 4, the previous equation reduces to equation (2). The minor dimension 
of noncircular nonplanar menisci should therefore be one -half of 6 e as given by equa- 
tion (1) for the same pressure P. 

In the case of an irregular hole S T e as given by equation (2) will generally apply. 
Temporarily, portions of the meniscus at the edge may be oblique with the critical minor 
dimension exceeding 8* . Consequently, the angle at the edge may be excessive and 
liquid escape will occur until the meniscus has retracted slightly into the capillary and 
acquired a critical minor dimension more nearly equal to or less than 6 T e . 

For a capillary radiator consisting of a panel of stacked parallel layers of uniformly 
spaced plates, equation (2) applies with the plate spacing 6 as the critical dimension. 
However, this type of capillary would have serious limitations unless internal edges 
were provided for liquid interface attachment. Suitable capillary structures would have 
further partitioning of the space between plates, or a cellular matrix included between 
the plates. Such a structure, in addition to being mechanically sounder, should enhance 
capillary retention by introducing additional surfaces and edges to which liquid may ad- 
here. A wick- or sponge -like structure, such as a foam metal matrix with intercon- 
necting pores, would be appropriate. This type of matrix, sandwiched between layers 
of sheet metal, should be homogeneous with characteristic "pore diameters" that lie 
within close bounds. Alternatively, the capillary medium may consist of parallel cir- 
cular tube segments in which case the tube diameter is related to surface tension 
and pressure through equation (1). 

System i nteraction factors . - The capillary phenomena upon which the radiator con- 
cept is based cannot be considered apart from the total system. Factors that influence 
capillary retention in a radiator loop under flow conditions become apparent by reference 
to the system depicted schematically in figure 5. The entire radiator loop, including all 
noncapillary volumes such as the pump, heat exchanger, and headers, is liquid -filled 
to capacity. In figure 5 a puncture is shown as a clean break in a cylindrical capillary 
tube that represents one of the numerous flow channels within the radiator panel. Be- 
cause the channels are interconnected, alternative flow passages invariably exist around 
each break as indicated. The radiator flow area may be taken as being very large com- 
pared to the area lost because of punctures. 
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Figure 5. - Schematic of all-liquid radiator loop with break in capillary channel. 


When a puncture occurs, a quantity of liquid is expelled along with a quantity of 
radiator material. If the lost liquid is not replenished, the resultant liquid interface 
will form within the capillary as in a sponge that is partially emptied by squeezing. 
Accordingly, menisci that form at the break site will be more or less displaced into the 
capillary channels. The menisci in figure 5 are both shown convex corresponding to an 
internal edge angle 9 greater than 90° and positive pressure. 

Given the aforementioned conditions, in order for liquid to escape, the menisci in 
figure 5 must move to the lip of the tube and the angle 9 must exceed a critical value 
greater than 90° at the lip to overcome the edge effect. Also, liquid must be withdrawn 
from elsewhere in the radiator since the amount of liquid is assumed to be constant. If 
a previous hole exists elsewhere, then the meniscus there must retract to supply the 
liquid. If the previous menisci have assumed curvatures that agree with the local pres- 
sures, no displacement can occur. If there are no previous holes, or if previous 
menisci are stable, then new menisci must form by a process that does not draw upon 
liquid from other parts of the system (unless additional liquid is supplied from a sep- 
arate reservoir. ) The capillary structure should be such that menisci can form at 
punctures sites with curvatures that agree with local pressures established prior to the 
puncture. 

Pressure rise in the radiator pump and absolute pressure levels in the radiator can 
be tailored to prevent potential leakage from punctures. Figure 6 presents a theoretical 
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(a) Wetting liquid, 0° <0 C <90°, 
A Pmax- r <sin e c + cos e c ). 


(b) Nonwetting liquid, 90° <0 C <180°, 

A Pmax- £<-«* 180°). 
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Figure 6. - Scheme for establishing absolute pressure and pressure rise across radiator pump. 


scheme for establishing the allowable pressure rise in the pump. In this example, 
hypothetical "pressure control openings" are located near the pump inlet and outlet. 

When the pump is operating, the inlet and outlet menisci are displaced as indicated. 
Provided the maximum allowable pressure rise AP is not exceeded, the menisci 
will be stationary. 

If the liquid wets (fig. 6(a)), then the maximum pressure at the outlet would be 
(4cr/5 e )[-cos (9 C + 90°)] and the minimum pressure at the inlet would be (4a/6 e )(-cos 0 c ). 
Recall that 9 Q + 90° equals 9 maximum. (The openings are cylindrical tubes with 
diameter 5 e and the external pressure is zero. ) This gives a maximum pressure rise 
of 


AP max = T (Sin + C0S °c ) (3) 

5 e 

Note that 1 < (sin 0„ + cos 9 J =£ 1. 414 for 0° < fl < 90°. 

If the liquid is nonwetting, then the inlet pressure control opening (fig. 6(b)) might, 
for example, be conical with a slope angle equal to 9 C so that the inlet equilibrium 
meniscus is flat and the corresponding liquid pressure is zero. At the outlet opening 
9 may not exceed 180° because the liquid pressure is maximum when 9 = 180°. For 
nonwetting liquids the maximum pressure rise may be at least 4a/5 e for all 9 Q . These 
cases are included to illustrate how the lowest pressure can always be made to assume 
values such that the allowable pressure rise AP max is at least 4a/ for all contact 
angles 0 . 

Equation (1) can be used to set the maximum pump or radiator pressure. Also, 
according to the examples given, -cos 0 in equation (1) may be set equal to 1 to get a 
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conservative estimate of AP m „ v where AP mov = P. Once the maximum pressure is 

IIlcLX XllcLX. 

specified by means of equation (1), a suitable minimum pressure can be obtained. 

Equation (3) exemplifies a case where AP can equal or exceed 4 c/6 a , but refer- 

max c 

ence 8 suggests that there might be practical limitations on stability when AP ex- 
ceeds 4a/6 e< 

The pressure profile of the radiator system established according to the previous 
hypothetical scheme should remain unaltered when the radiator is punctured. As ex- 
plained previously, this can be assured if the capillary matrix has a substructure that 
provides a suitable array of surfaces and edges. 

Liquid drainage at puncture sites in the radiator would tend to occur if noncapillary 
components such as the pump, headers, or heat exchanger were partially voided and 
subjected to a sufficiently high pressure. The absence of such voids is therefore a pre- 
requisite. When filled to capacity, systems similar to that shown schematically in fig- 
ure 5 will retain liquid because of cohesive forces in the liquid in accordance with the 
previous argument. 

If capillary pumping is used, a vapor void is an integral part of the pump. Capillary 
pumping depends in part on vaporization and condensation for fluid transfer. In this case 
vapor -liquid interfaces arise in the capillary matrix within the pump. A radiator loop 
with a low pressure vapor void of this type is capable of holding liquid because of inher- 
ent capillary retention within the pump as well as at puncture sites. 

The effects of inertial pressures that arise from vehicle accelerations and pulse 
pressures that accompany high velocity impacts would engender liquid escape if large 
pressure fluctuations result. Although inertial pressures can be taken into account by 
design, it is difficult to assess the effect of pulse pressures without experimental data. 


Heat Transfer and Fluid Dynamic Considerations 

Parametric equations relating heat transfer, hydraulics, fluid properties, and radi- 
ator geometry are derived for a hypothetical capillary radiator in the following sections. 
The purpose of these sections is to evolve analytic expressions by means of which the 
concept may be evaluated. The evaluation is based on the possibility of selecting a 
capillary dimension consistent with liquid retention at meteoroid impact sites under the 
constraints imposed by heat transfer requirements. 

Identification of conditional parameters. - In the ensuing discussion the radiator 
consists of a number rip of identical thin, flat heat radiating panels. The panels are 
joined to a central header system and are so disposed with respect to one another as to 
form the dihedral angle o> as depicted in figure 7. The total radiating area equals 
2LWn where L is the length and W is the width of each panel. Edge radiation is 

it 


10 




Figure 7. - Capillary radiator with eight radiating surfaces. 


ignored because of the small area involved. Given Q, the amount of heat to be rejected, 
the required radiator area is 


2LWru = 9 (4) 

P q 

The radiative heat flux q is a function of the surface emittance, temperature, and 
dihedral angle between panels (the space environment being regarded as a perfect heat 
sink); that is, 


* = V ct B T s (5) 

In the present case the emittance e and T g , the radiator surface temperature, are 
assumed to be constant over the radiator surface. The M configuration factor" f as- 
sumes values that vary with the dihedral angle oj. An equation relating f to a> is 
derived in appendix B. 

The flow in the panels is shown in figure 8. The radiator panels consist of alternat- 
ing layers of fluid and metal as depicted in figure 9. Each fluid layer is further parti- 
tioned into parallel capillary channels extending across the length of the panels. In each 
panel the capillary channels accept fluid from a common supply header and empty into a 
common return header (fig. 8). Header arrangement and size do not enter explicitly into 
the analysis that follows. The channels are square tubes and have a "nominal diameter" 
equal to 8. 
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Figure 8. - Schematic of flow scheme in radiator panels. 
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Figure 9. - Representative capillary medium contained within radiator panels. 



There are layers of channels which, combined with (n^ + 1) alternate layers of 
surface and partition material, make up the radiator panel thickness in accordance with 
the following equation: 


t r - n^5 + (n^ - l)t c + 2tp 


( 6 ) 


Note that the radiator surface plates and partitions are not necessarily identical in thick- 
ness, t ^ t . The quantity n.5 is termed total fluid layer thickness t 7 ; that is, 

P C L L 

h = V (7) 

There are n £ parallel channels in each layer so that the length L is given by the 
equation 


L = n c 5 + (n c + l)t c s n c (5 + t c ) (8) 

The approximate expression is used herein since n c is considerably greater than one. 

Heat transfer and flow relations . - Conductive heat transfer within the radiator 
delivers heat to the radiator surface at a rate equal to q as is indicated by the following 
equation: 


q = 


2k(T i 


T ) 
s' 


(9) 


In equation (9) k is an average value representing the combined thermal conductivities 
of the containment structure and fluid, and is the maximum internal temperature. 
The mass flow rate in each panel M is 


M = 


Q 


v AT 


( 10 ) 


where AT is the temperature difference between the radiator inlet and outlet. The 
mass flow rate is also given by equation (11) 


M 


P vn c n l 


(U) 


The flow velocity v is obtained by the Hagen-Poisseuille equation for laminar flow in 
circular tubes which is applicable in this case; this equation is 
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AP = 1 27 ? Wv (12) 

6 2 

Arbitrarily setting AP equal to Ip and solving for v yield equation (13). (This pro- 
cedure assumes that a pressure of ~P is sufficient to overcome the head loss in the 
return header. ) 


v = 


_P6__ 

64r?W 


(13) 


General parametric equation . - Setting 5 equal to 6 e , using equation (1) for this 
case, combining equations (1), (4), (8), (10), (11), and (13), and solving for L in terms 
of 6 yield the following parametric equation: 


L = 



(14) 


where 


= 8Q 2 

AT(-cos #)qn 2 

r 


(15) 


and 





(16) 


A useful form of equation (14) is obtained by combining with equation (7) and rear- 
ranging slightly to obtain 


l 



(17) 
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Equation (14) incorporates all the pertinent quantities associated with the radiator save 
those involving conductive heat transfer as expressed by equation (9). The subscript Tj 
appended to the bracket in equation (16) indicates that the values assigned to 77 , c, p, 
and a must be those that are appropriate to the radiator internal temperature. 

In the case of the header arrangement shown in figure 7 the quantity L, represent- 
ing both header and panel length, imposes a constraint upon 6. An excessive header 
length, for example, necessitates high pressure drops. This implies, in turn, high capil- 
lary pressures that may exceed that required by equation (1) for a given capillary diam- 
eter. A significant connection between header length and possible inertial pressures is 
discussed in appendix C. 

In equation (15) the thermal power of the radiator Q is fixed by the mission and/or 
system requirements. Radiator temperatures T g , Tj, and AT are selected to assure 
that the properties of the heat exchange liquid assume values compatible with the radi- 
ator requirements; in particular, Tj is selected to produce a low vaporization rate. 

The radiative heat flux q depends upon the nature of the radiator surface selected, the 
surface temperature T_, and the dihedral angle w, which is determined, in turn, by the 
number of panels, n The quantity -cos 9 may be arbitrarily set equal to a value rang- 
ing from 1/2 to 1 corresponding to 6 < 120° to 180°. The parameter is thus a con- 
stant that is predetermined by the operating conditions and physical and geometric prop- 
erties of the radiator. 

The parameter \p, which is a composite fluid property, is useful for rating liquids 
contemplated for use in the capillary radiator. Liquids having low \p are preferred 
because, as evident from equation (14), the result is smaller header lengths L in fig- 
ures 7 and 8, and hence less header pressure drop. Note that \p has the dimensions of 
a reciprocal heat transfer coefficient. Fluids with lower values of xp have better heat 
transfer characteristics. 

Parametric analysis procedure. - Once the parameters 4> and ip have been as- 
signed values, equation (14) may be used to investigate the range of possible practical 
capillary diameters. An upper limit is imposed on 6 through equation (1) by limits on 
the minimum operating pressure P, that can be attained. A second upper limit is im- 
posed as a consequence of inertial pressures that result from vehicle accelerations as is 
explained in appendix C. A third upper limit on 6 is imposed by adoption of the Reynolds 
criterion as follows. 

The applicability of the Hagen-Poisseuille equation used in obtaining equation (14) 
depends on the critical Reynolds number for capillary channels. According to refer- 
ence 9, experiments have indicated that Darcy’s Law and, consequently, the Hagen- 
Poisseuille equation applies to capillary channels for Reynolds numbers less than 
roughly 75. This criterion is not well established, and may be several orders lower 
for capillary channels formed in some types of porous media. Herein the critical 
Reynolds number is taken as 100; that is, 
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Re = < 100 

V 


( 18 ) 


Combining with equation (13) and substituting P from equation (1) yields 

Re = p qj . -cos 6) (^\ < 100 (19) 

16 r, 2 V W / 

The mi nimum value of 6 is determined by considering the effects of reducing 6 
upon thermal conduction and radiator weight. Examination of equation (17) reveals that 
as 6 decreases t^ increases for fixed values of <&, \js, and t c /6. Hence, as is evi- 
dent from equations (6) and (7), decreasing 5 results in increasing t r . An increase 
in t is necessarily accompanied by a corresponding increase in the internal tempera- 
ture T- and a lesser increase in the surface temperature T (see eqs. (5) and (9)). 

1 

However, T. cannot be increased indefinitely, for the internal temperature must be in 
a range that does not result in excessive vapor pressure nor extremely low surface ten- 
sion in the radiator fluid. Hence, a lower bound on 6 is set by the overall panel thick- 
ness and thermal conduction. 

Radiator specific weight also influences the selection of a minimum 5. Note that 
partitioning material must be added for each diminution of 6; that is, smaller capil- 
laries entail higher containment weight. The consequence of reducing 6 is an increase 
in n^ and t^ (eq. (17)) and a corresponding increase in radiator specific weight with 
respect to generated power w. This is apparent from examination of equations (21) and 
(22) in the next paragraph. 

An expression for specific weight is obtained by referring to figures 8 and 9 and 
taking into consideration both the weights of the radiator panels and contained liquid: 

w = w £ + Wj (20) 

where 



and 
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(22) 



In equation (21) the material densities of the surface plates and capillary partitions 
represented by thicknesses t and t , respectively, are assumed equal. Equation (21) 
is obtained by using the approximate form of equation (8), given that n Q » 1. 

As 6 assumes smaller values, either t may be constant or the ratio (t c /6) may 
be constant. In each case, all the other quantities in equations (21) and (22) being fixed, 
both w c and w^ increase with n^6(=t^) which, in turn, increases with decreasing 6 
as stated previously (see eqs. (7) and (17)). 

The selection of heat exchange liquid, mean radiator temperature, and other quan- 
tities may be predicated on the margin between the minimum and maximum 6. The co- 
incidence or overlap of maximum and minimum values of 6 would indicate impractica- 
bility under the conditions selected. Because of the monotonic logarithmic functions 
involved, as in equation (14), the maxima and minima referred to previously are not 
analytical bounds, but are dependent upon the relative limits assigned to factors such as 
specific weight, pressure, structure dimensions, and so forth. In general, it is de- 
sirable to make 5 as small as possible consistent with the limitations mentioned before. 


Liquid Loss Considerations 

Meteoroid impacts on a capillary radiator panel will result in pressure pulses and 
punctures. There is insufficient relevant information concerning the intensity and trans- 
mittal of pulse pressures in liquid-filled capillary media for estimating liquid escape 
caused by impact pressures. However, liquid loss by displacement or temporary efflux 
preceding stabilization of the meniscus, and vaporization can be estimated by consider- 
ing the frequency and geometry of punctures. Equations for computing the decrement of 
radiator surface and flow channel areas and liquid loss by displacement and vaporization 
are presented in this section in correspondence with the radiator configuration described 
in the section Identification of co nditional parameters . 

Meteoroid env ironment . - To determine the effects of the meteoroid environment, it 
is necessary to select a meteoroid "mass -frequency" distribution relation that best rep- 
resents the situation during the mission. Mass -frequency distribution data lie in two 
general categories, one applying to near -Earth and the other to "deep space" or inter- 
planetary regions. Accurate data for deep space are scarce, and curves representing 
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frequency distributions in the vicinity of the Earth disagree somewhat. However, an 
analytical expression needs to be adopted in order to estimate the probable rate of mete- 
oroid encounter and extent of damage that may result. A widely accepted meteoroid 
mass -frequency curve is that of reference 10 known as the Whipple -1961 curve. Al- 
though a later revision, Whipple-1963 (ref. 11), indicates a slightly lower meteoroid 
mass frequency, the earlier version as expressed in the following equation is adopted 
here: 


F = 1. 246X10 -15 m -1 " 1 (23) 

Equation (23) represents the cumulate frequency of all meteoroids having a mass equal 
to and greater than m. 

An expression for the total mass encountered by a unit area per unit time for the 
mass range m^ to mg is obtained as follows: Let dF m be the total mass flux in the 
range dm about m. Then. 


dF m = F dm (24) 

Integrating between limits m^ and mg yields 

F m = / ^ F dm, mg > nij (25) 

•'m^ 

Equation (25) implies that particles with masses less than m^ are inconsequential while 
those with masses greater than mg are too infrequent to be considered in lieu of the 
radiator area and duration of the mission contemplated. 

It should be noted that the aforementioned differs from the approach usually followed 
in determining anticipated meteoroid encounters (e. g. , ref. 1). Conventionally, a 
Poisson distribution is assumed in order to arrive at an expression for the probability 
of a vulnerable component being hit by a mass greater than m. Herein it is assumed 
that the number of encounters for any mass range equals the time -area average of the 
integrated flux F, that is, a probability factor of one is assumed. This approach places 
a greater burden on the capillary radiator concept. The survivability of the capillary 
radiator, unlike the fin -tube radiator, is not rated in terms of the probability of a punc- 
ture and consequent total disablement of a vital component. For a puncture in a capillary 
panel is not equivalent to total loss of that component. 

Puncture criter ia and loss equation s. - All ’’complete" punctures are idealized as 
cylindrical holes passing through the radiator panel (fig. 10). A puncture of one surface 
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Figure 10. - Cross section of idealized circular meteoroid puncture site. 


plate on either side of a panel is considered a complete puncture. Partial penetrations 
of a surface plate are ignored because the previous criterion for complete punctures is 
deemed broad enough to include all significant penetrations. Because of the results of 
the experiments mentioned in appendix D, the diameter of the hole is considered to be 
invariably twice the "diameter" of the impinging meteoroid irrespective of its shape, 
mass, or velocity. The meteoroid is assumed spherical, and its diameter d m is there- 
fore determined by its mass and density; that is, 


m 



(26) 


and the radiator area decrement for each puncture is 


a = 



(27) 


The differential fractional surface area decrement per unit time dA g , resulting 
from puncture by masses in the range dm about m, is 


^s 



(28) 


Using equations (26) and (27) to solve for a as a function of m, substituting for a, and 
F from equation (23) in equation (28), and integrating yield 


A 


s 


(-)3. 637xl0“ 14 p m “ 2 / 3 m“ 1 / 3 


m 2 

m l 


(29) 
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(The exponent of m in eq. (23) is equated to 1 for simplicity. ) The fraction of radiator 
surface area lost during the mission is 


A ts = tA s 


(30) 


where r is the mission duration. 

With each puncture the flow channel area diminishes. Because the capillary channels 
are interconnected, the flow is ’’two-dimensional”; therefore, the fractional loss of 
hydraulic area A rh is directly proportional to half the fractional loss of radiator area 
A _ if a random puncture pattern exists. The factor one -half appears because A 0 in 
equation (30) represents the two surfaces between which the capillary medium is in- 
cluded. 

Similarly, equation (30) also provides an estimate of the fraction of liquid lost by 
displacement when a puncture occurs. If the diameter of the slug of displaced liquid 
equals the puncture diameter, then D^, the fractional liquid inventory loss by displace- 
ment, is given by the following equation: 


A 


D . =■ 
rZ 


TS 


2x 


(31) 


In equation (31), x is the factor by which the liquid content of the entire radiator system 
exceeds the content of the panels. 

In figure 10, the area of the hole in the radiator surface and not the exposed liquid 
area is arbitrarily considered the ’’effective” area for evaporation. The actual ’’effec- 
tive” area depends partially on t and the depth -to -diameter ratio of the liquid inter- 
face. However, the approach used fixes the area for efflux even if the interface has re- 
tracted into the capillary channels. Accordingly, the vapor efflux area for a puncture is 
permanently established by equation (27). The increment in vapor efflux area caused by 
masses in the range dm about m equals dA„ (eq. (28)). Assuming a linear increase 

o 

in the vapor efflux area with time, the average area for vaporization accrued during the 
mission is 


1 QA S 

A = - n LWA = — ^ 
rZ 2 P TS 4q 


(32) 


The rate of evaporation is given by Langmuir’s equation (ref. 12) as 


x-1 


/M„ 


G = 5.833X10 A v + 

mm 1 / ^ 

A i 


(33) 
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The initial weight of liquid in the radiator, headers, and reservoir is, approximately, 



The fraction of liquid inventory lost during the mission is 


Ti h 


or, substituting from equations (14), (32), (33), and (34), 


( 34 ) 


(35) 


I j = (-)l. 0605X10 -14 





■ m. 


'm. 


(36) 


Note that equation (36) provides a criterion for selecting the heat exchange liquid and 
radiator temperatures, Tj and also T g , if 1^ is preselected. 

A further consideration is the relation between the surface plate thickness t 
and the depth to which an impinging meteoroid penetrates the surface plate. The rela- 
tion depends on the ratio between depth of penetration p and meteoroid diameter d , 
that is, p/d m . If p/d , the meteoroid velocity, mass m, and density p m are speci- 
fied, it is then possible to estimate the surface plate thickness in which only partial pen- 
etrations are likely to occur for meteoroid masses less than a given value. A suitable 
criterion for selecting t should require that for masses less than a minimum value, 

Jb' 

for example m^, the effect is equivalent to no penetration at all. Surface plates thinner 
than this would result in higher rates of liquid exposure because of the higher frequen- 
cies of encounter with smaller masses predicted by equations that describe the mete- 
oroid environment (e. g. , eq. (23)). Therefore, the selection of tp depends on the 
allowable penetration threshold on one hand and its effect on the radiator weight on the 
other (see eq. (21)). 

Establishment of the criterion for tp mentioned previously is difficult because p 
and hence p/d m varies with the meteoroid mass and velocity and because meteoroid 
densities and shapes vary widely (ref. 13). However, the following approach can be 
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used: set t p equal to p and assume a value for p/d m corresponding to the mean 
probable velocity. Consequently, p and hence tp may be regarded as a function of the 
meteoroid mass through equation (26). Using equations (30) and (36), t (=p) may be 
related through m to the fractional area and liquid inventory losses. Selection of the 
surface plate thickness will be based on that value of t for which the maximum accept- 

r 

able area or liquid loss occurs. 

DISCUSSION OF APPLICATIONS AND PROBLEMS 
Illustrative Mission, Configuration, and Conditions 

The potential feasibility of the capillary radiator concept may be assessed in terms 
of the latitude available for selecting such crucial quantities as the capillary diameter, 
panel dimensions, temperature, hydraulic pressure, and other pertinent dimensions. If 
a compatible set of the aforementioned quantities exists, it is then necessary to deter- 
mine whether they lead to acceptably low levels of specific weight, area attrition, and 
fluid loss. An assessment of potential feasibility in terms of the criteria derived in the 
previous sections is presented hereinafter in connection with a specific radiator applica- 
tion. 

The illustrative example developed in this section concerns a prospective applica- 
tion of the capillary concept which affords a means of comparing the specific weights of 
capillary and fin-tube radiators. The total radiated heat rejection is arbitrarily set at 
approximately one -third megawatt. This power level is in the general range contem- 
plated for power system cooling in manned interplanetary vehicles and lunar base power- 
plants (refs. 14 and 15). Smaller versions of capillary radiators may have applicability 
to subsystem or component cooling in space vehicles. Examination of equations (14), 

(15), (21), and (22) will show that for lower power levels the vulnerable area and specific 
weight are also lower. 

The radiator configuration for the following example consists of four orthogonal 
panels joined to a common header as illustrated in figure 7 (p. 11). Fluid flow within the 
panels is perpendicular to the header axis as illustrated in figure 8 (p. 12). The capillary 
substructure is the same as that described previously and illustrated in figure 9 (p. 12). 

The following conditions are prescribed: 


Duration of mission, r, sec (days) 8. 64x10^ (1000) 

Power cycle efficiency, e 0. 333 

Radiator power (rejected heat), Q, W 0. 3xl0 6 

Radiator surface temperature, T s , °K 589 

Internal fluid temperature, Tj, °K 655 

Radiator power flux, q, W/m^ 4. 3xl0 3 

Heat exchange fluid Liquid lithium 

Material of construction Titanium 
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Radiator power rejection is based on a cycle efficiency of 33 percent. To evaluate 

equation (15), the radiator inlet to outlet temperature difference AT is 55. 5° K, and for 

simplicity it is assumed that all heat is radiated to space at a uniform radiator surface 

temperature T„ of 589° K. The temperature level of the radiator is arbitrarily se- 
s 3 

lected to yield a low vaporization rate at puncture sites. The heat flux of 4. 3x10 watts 

per square meter is based on a radiator surface emissivity e of approximately 0. 9 and 

a configuration factor f of 0. 8 (see appendix B). The value of 9 in equation (15) is 

taken as 120 degrees so that -cos 9 = 1/2 ^ 1, the pressure drop criterion cited in the 

section System interactio n factors . 

Lithium is selected as the heat transfer liquid primarily because of its low vapor 
pressure in the temperature range contemplated in the example. Titanium is selected 
as the containment material because of its compatibility with lithium and its low density. 
However, neither the fluid nor containment material need be specified. The crucial 
properties such as density, vapor pressure, and the quantity xfy can be independently 
specified or treated parametrically. Further analysis might then point to suitable sub- 
stances. 

By selecting the particular radiator configuration, materials, and conditions men- 
tioned previously, it is possible to concentrate upon three factors that are significant in 
evaluating the capillary radiator’s feasibility: capillary diameter, radiator surface at- 
trition, and liquid loss. 

Bounds on capilla ry dimension . - The selection of the capillary diameter is based on 
equation (1) in this case and evaluated in accordance with the discussion in the section 
Parametric analysis p rocedure . Substitution of the conditional values for <3> and i p in 
equation (14) yields 


L = 


1.408X10" 2 



(37) 


A plot of L against 5 appears in figure 11 where t c /6 is set equal to 10"^ and n^ 
assumes the values 1, 2, 4, 9, 25, and 81. In this example i f/ and a are assigned 
values corresponding to T^ = 655° K for lithium to obtain the expression for L in 
equation (37). 

According to the discussion in the section Parametric analysis pro cedure , upper 

bounds on 6 are imposed by the Reynolds number, hydraulic pressure requirements, 

and inertial forces. In the present case these factors collectively dictate an upper bound 

-3 -2 

on 5 of about 1. 5x10 to 10 meter. The calculations for this and the following are 
discussed in appendix E and the results are summarized in table I. 
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Figure 11. - Covariation of radiator length with capillary dimension 
and number of capillary layers according to equation (14) for 
3xl(r-watt radiator. 


TABLE I. - EFFECT OF CAPILLARY 


DIMENSION ON SPECIFIC WEIGHT 3 - 


Capillary 

diameter, 

8, 

m 

Number of 
layers, 

Specific 
weight of 
panel alone, 

w c> 

kg/kW 

Specific 
weight wit 
liquid, 

w, 

kg/kW 

1. 5X10" 3 

1 

0. 67 

0.82 

1.1 

2 

.84 

1. 1 

.73 

4 

1.0 

1.3 

.49 

9 

1.4 

1.8 

.29 

25 

2.0 

2.8 

. 16 

8! 

3.2 

4. 5 


a Based on equations (21) and (22) and lithium- 


filled titanium panels with radiator panel 
and header length L = 10 m, radiator sur- 
face temperature T g = 655° K, thickness 
of surface plate t = 2. 54x10"^ m, t /6 = 
0. 1, power conversion system cycle effi- 
ciency e = 0. 333, and radiator heat flux 
q = 4.31X10 3 W/m 2 . 
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The lower bounds on 6 are imposed by heat transfer and specific weight limitations. 
Given the dimensions and materials of this example, it is found that the more significant 
lower bound is imposed because of the increase in specific weight w that accompanies 
the reduction of 5. For instance, 5 has a lower bound of roughly 1. 5x10”^ meter if a 
maximum specific weight of about 5 kilograms per kilowatt is prescribed. 

To satisfy the criteria mentioned in the previous paragraphs, 6 should lie in the 
-4 _3 

range between 1. 6x10 and 1. 5x10 meter, approximately. These figures are in the 
range in 6 from 1. 5x10”^ to 1. 5Xl0 - ^ meter associated with n^ from 81 to 1, respec- 
tively, for the arbitrary selected length of L = 10 meters (see fig. 11). 

In the case presented herein it is inferred that the capillary radiator concept is 
feasible in terms of the heat transfer and fluid dynamic factors corresponding to the 
bounds on 5. It is noteworthy that the specific weights of a 1/3 -megawatt capillary radi- 
ator having between 1 and 81 liquid layers are in the same range as those quoted in ref- 
erence 4 and summarized in table II for a 1/2-megawatt conventional fin -tube radiator. 
The basis of comparison in this instance is solely the specific weight of the radiator 
panels exclusive of headers, pumps, and other appurtenances. Note that the specific 
weights in table II are the minimum values quoted in reference 4 for each material. 

Surface area and liquid losses . - Rough estimates of radiator surface area attrition 
and liquid loss may be made by means of equations (30), (31), and (36). Using the con- 

TABLE n. - MINIMUM SPECIFIC WEIGHT FOR 


1/2 -MEGAWATT FIN-TUBE RADIATOR a 


Armor materials for 
tube protection 

Specific weight, 
kg/kW 
(b) 

Columbium alloy (Cb-lZr) 

5 

Nickel alloy (Inconel 718) 

3. 7 

Cobalt alloy (L-605) 

3. 6 

Steel (SS 316) 

3.3 

Molybdenum alloy (TZM) 

2.8 

Vanadium 

2.6 

Graphite 0 

1. 3 

Beryllium 0 

1.0 


a From ref. 4, specific weights of empty 
tubes. 

^Minimum specific weight occurs for all 
materials at radiator temperatures 
between 920° and 1030° K in fin -tube 
radiator. 

c Includes thin-waUed inner liner tube of 
columbium alloy (Cb-lZr). 


25 



TABLE IB. - RADIATOR AREA LOSSES AND LIQUID 
LOSSES FROM METEOROID PUNCTURES 3 ' 


Meteoroid 

Surface 

Fractional 

Fractional 

mass, 

plate 

radiator 

liquid 

mi, 

thickness. 

area loss, 

inventory 

kg 

t b 

V 

A rs or D rl C 

loss, 


m 


1 7 d 
Tl 

10" 4 

2. 5X10" 1 

1. 7x10' 7 

0. 7xl0' 6 

10' 6 

2. SxlO' 2 

3. 0xl0” 6 

1.4X10" 5 

10" 9 

2. 5X10' 3 

3. IxlO' 5 

1. 5xl0' 4 

10- 12 

2. 5xl0' 4 

3. 1X10' 4 

1.5X10' 3 

10' 15 

2. 5X10~ 5 

3. lxlCf 3 

1. 5X10' 2 

10- 18 

2. 5xl0~ 6 

3. IxlO' 2 

1. 5X10' 1 


a Estimates based on eqs. (26), (30), and (36) 


for 1000 -day mission. 

D tp equated to p (a function of m^). 

c 1 

Multiply by | to obtain fraction of liquid lost 
by displacement D 7 . 

a Based on characteristic capillary dimension 
6=1. 5xl0 -3 m. 

ditions given at the beginning of this section, it was found in appendix E that a 1000-day 

mission should entail acceptably small losses of area and heat exchange liquid. As in- 

_ 4 

dicated in table in, for example, with the surface plate thickness tp = 2. 5x10 meter, 
the estimated surface area attrition is 0. 03 percent, and liquid loss by evaporation is 
0. 2 percent, approximately. 

Also note that the previous liquid loss figures support the radiator temperature 
range selected, that is, 589° to 655° K. In fact, for lithium any temperature from ap- 
proximately 589° to 868° K results in a low loss rate since the corresponding values 
of G, a, and if/ do not change significantly for this temperature range. An increase in 
the radiator temperature permits a corresponding increase in q, the radiative flux, and 
a corresponding expansion of the bounds of 5 and associated quantities through equa- 
tions (14) and (15). 

The figures quoted for the bounds on 6 and for area and liquid losses indicate 
feasibility at least under the selected conditions described in this example. The example 
given also indicates that a capillary radiator is potentially competitive with conventional 
radiators up to and including the megawatt range. 



Problem Areas and Prospective Advantages 


A definitive assessment of the feasibility and applicability of the capillary radiator 
concept must be deferred until the experimental verification of certain crucial phenomena 
has been accomplished. For example, studies of surface tension and vaporization phe- 
nomena in capillary media under high temperature, vacuum conditions are needed to 
support the values assumed herein. For the same reason, specific and accurate infor- 
mation concerning heat transfer, flow, pressure drop, and temperature distribution in 
porous media of the types suitable for a capillary radiator is needed. Indeed, the pro- 
duction and fabrication of capillary media that can be effectively incorporated into radi- 
ator panels are potentially problematic and require developmental effort. 

Some aspects of capillary retention under dynamic flow conditions may require scale 
model studies for further verification of workability. Also, the possible presence and 
effects of noncondensable gases on capillary retention should be assessed. 

A vital area for further inquiry concerns the effects of meteoroid impacts on liquid - 
filled capillary media. Experimental information is needed to define the characteristics 
of impact damage in various capillary structures and to determine the effects of intense 
pressure pulses on the structures and contained fluid. For example, under certain con- 
ditions, highly localized, short duration pressure pulses may be exerted which are of the 
order of billions of newtons per square meter and last roughly 100 microseconds 
(ref. 16). The evolution of suitable capillary substructures and modes of fabrication 
will be necessary to circumvent the possible deleterious effects of impact shocks. The 
aforementioned experiments and developmental efforts should be tailored for testing 
under conditions duplicating the zero - or low -gravity environment of an actual space 
application. 

Because the heat exchange liquid should have high surface tension and low vapor 
pressure, the choice of fluids for high temperature operation may be somewhat re- 
stricted. For example, only liquid metals appear to possess suitable properties in the 
high temperature, kilowatt to megawatt range. In this case, the problem of containment 
material compatibility also arises. An investigation of corrosion and dissolution phe- 
nomena in porous media would be required in order to establish evidence of long term 
integrity of capillary media. In addition, means of inducing wetting where needed while 
inhibiting unwanted spreading of the liquid over the radiator surface should be developed 
to optimize the surface tension phenomena upon which the retention concept is predicated. 

Incorporation of the capillary radiator concept into a viable system poses a number 
of engineering problems. The evolution and development of suitable modes of pumping, 
fluid distribution; and pressure, temperature, and inventory control may be problematic. 
These factors might dictate the ultimate forms, applications, and limitations of the cap- 
illary radiator concept. 
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Among the potential advantages is the possible reduction of radiator weight with re- 
spect to power. Excluding the pump, headers, and supports structure, the radiator 
specific weight might be reduced substantially below values usually quoted for armored 
fin -tube radiators. Another potential advantage is afforded by the low operating tem- 
perature which may be as much as 650° K lower than that contemplated for an equivalent 
fin -tube radiator (refs. 1 to 3). The temperature advantage accrues because a consider- 
ably larger, prime heat radiating surface area is possible. This, in turn, is a conse- 
quence of the fact that it is necessary neither to confine the heat exchange liquid in a 
small number of widely spaced flow channels nor to depend crucially upon conductive 
heat transfer to the radiating surface. The low radiator temperature may be utilized 
either to increase cycle efficiency or to lower the overall temperature level of the power 
system. Also the overall size or area of the radiator is not as critical as in a conven- 
tional fin -tube radiator because the area vulnerable to punctures need not be reduced to 
avoid punctures. 

In addition to the potential advantages mentioned previously, it may be possible to 
eliminate redundant radiator segments that are contemplated for fin -tube radiators 
(ref. 6). Also, the thin panels or fins in a capillary radiator may offer more configura- 
tional flexibility for launch folding and deployment. 

It should be noted that prospective applications of this concept are not limited to heat 
rejection in large powerplants, but may include the cooling of smaller components and 
subsystems. 


CONCLUDING REMARKS 

The capillary radiator concept described herein is suggested for use in space heat 
rejection systems. If the requisite conditions for capillary retention of liquid described 
herein are met, it appears possible to construct capillary radiators that satisfactorily 
survive the effects of meteoroid impact damage. 

Although liquid retention at puncture sites by surface tension effects is theoretically 
possible by using appropriate capillary media, more studies are needed to confirm the 
viability of a liquid -filled capillary radiator under the full range of conditions prevailing 
in an actual application. Such effort seems merited because a capillary radiator may 
offer advantages that are suitable to space power systems. The chief potential relative 
advantages are (1) lower system temperature or higher theoretical powerplant cycle 
efficiency, (2) adaptability to various power systems and configurational flexibility, 

(3) lower panel weight with respect to electric power produced, and (4) less restriction 
on the radiator size. 
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Before the potentialities of the capillary concept can be realized, however, a num- 
ber of technical and developmental problems will have to be resolved. In particular, the 
effects of hypervelocity impacts on liquid-filled capillary media should be investigated 
and means for producing suitable capillary structures should be developed. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 25, 1967, 

120-27-02-08-22. 
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APPENDIX A 


SYMBOLS 


A g fractional area loss per unit time 
(eq. (29)), sec -1 

fraction of hydraulic cross section 
lost 

A,^ mean equivalent liquid area exposed 
(eq. (32)), m 2 

A tS fraction of radiator area lost during 
mission (eq. (30)) 

a radiator area decrement, m 

c specific heat of radiator liquid, 

J/(kg)(K) 

fraction of liquid lost by displace- 
ment (eq. (31)) 

d puncture diameter, m 

d m meteoroid diameter, m 

e power conversion system cycle 

efficiency 

F cumulate meteoroid frequency, 

(m -2 )(sec - ^) 

F__ integrated mass flux (eq. (25)), 

kg/(ni )(sec) 

f configuration factor (eq. (38)) 

o 

G evaporation rate, kg/(m )(sec) 

g acceleration, m/sec 2 

g Q acceleration due to gravity, 

9. 80 m/sec 2 

1^ liquid inventory, kg 

I . fraction of inventory lost 

(eq. (36)) 


i ff quantity proportional to total heat 
flux when <jj = 180° 

i quantity proportional to heat flux 

for dihedral angle o> (eq. (39)) 

k thermal conductivity, W/(cm)(K) 

L radiator panel and header 

length, m 

di differential length of meniscus 

(fig. 4), m 

M mass flow rate per panel, kg/sec 
M gram molecular weight, g 

o 

m meteoroid mass, kg 

m^ minimum mass (integration 

limit), kg 

mg maximum mass (integration 

limit), kg 

n c number of capillary channels per 

layer 

n^ number of liquid layers 

rip number of radiator panels 

2 

P radiator inlet pressure, N/m 

2 

AP radiator pressure drop, N/m 

o 

P inertial pressure (eq. (43)), N/m 

o 

P mm va P or pressure at |(Tj + T g ), mm Hg 

p depth of penetration, m 

Q thermal power rejected by radi- 

ator, W 

o 

q radiator heat flux, W/m 
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Re capillary Reynolds number 
(eq. (19)) 

AT inlet to outlet temperature differ- 
ence in radiator, K 

Tj maximum radiator internal tem- 
perature, K 

T radiator surface temperature, K 

s 

t £ thickness of capillary or layer 
partition, m 

t^ total fluid cross section thick- 

ness, m 

tp thickness of surface plate, m 

t overall thickness of radiator 

panel, m 

v liquid velocity in radiator panel, 

m/sec 

W width of radiator panel, m 

w total specific weight of radi- 

ator panel (eq. (20)), kg/W 

w c specific weight of capillary panel 
(eq. (21)), kg/W 

w^ specific weight of heat exchange 

fluid (eq. (22)), kg/W 

x inventory factor (x > 1) (eq. (31)) 

y variable length (eq. (40)), arbi- 

trary units 

/3 variable angle (eq. (39)) 

6 characteristic capillary dimen- 

sion, m 

6 g equilibrium capillary dimension for 
tubular capillary (eq. (1)), m 


6g equilibrium capillary dimension for 
inertial pressure (eq. (44)), m 

5’ local capillary dimension 

(fig. 4), m 

6^ equilibrium capillary dimension 
for parallel plate capillary 
(eq. (2)), m 

6” arc length of differential meniscus 
(fig. 4), m 

e therma.1 emittance 

2 

rj viscosity, (N)(sec)/m 

9 0 c < e < 180°, angle with respect to 

container surface at edge, deg 
(fig. 2(b)) 

6' liquid surface angle of contact with 
respect to contiguous element of 
liquid, deg 

6 c equilibrium contact angle, deg 

O 

p liquid density, kg/m 

0 

p c channel material density, kg/m 

q 

p m meteoroid density, kg/m 

a surface tension, N/M 

<Xg Stefan-Boltzmann constant, 
5.6697X10" 8 W/(m 2 )(K 4 ) 

t mission duration, sec 

parametric constant (eq. (15)) 

cp variable angle (eq. (39)) 
fluid parameter (eq. (16)) 

co dihedral angle between radiator 

panels, deg 
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APPENDIX B 


CONFIGURATION FACTOR 


In a radiator consisting of panels that are arranged to form dihedral angles as in the 
case illustrated in figure 7, a portion of the radiated heat is intercepted by adjacent 
panels. The amount that is intercepted depends upon the dihedral angle co. The fraction 
of the total heat radiated that escapes to space (regarded as a perfect heat sink) is a 
function of the temperature, emissivity, reflectance, and absorptivity integrated over 
the surface of the panel. For simplicity, the configuration factor derived hereinafter is 
based upon the following conditions: 

(1) The radiation intensity in any direction varies as the cosine of the angle to the 
normal of the surface, Lambert’s cosine principle. 

(2) Temperature and emissivity are constant over the radiator surface. 

(3) The panels are infinitely long and equal in width. 

(4) Reflected radiation is nil, and all intercepted reemitted radiation is nil. 

The configuration factor f is defined by equation (38) 


f 

U) 



(38) 


where i is proportional to the total heat flux emitted by a panel that forms a dihedral 
angle co with respect to an identical adjoining panel and i is proportional to the total 
heat flux when ou equals 180° (see fig. 12). The defining equation for i is 

/ I /• n/2 rl T 

/ / sin fi dy d cp d/3 (39) 

v J 0 pj 

where 

ff-. = arctan / '- in — cos ? ) (40) 

x \ y - cos co / 

Because the width of each panel W is identical, the limit on y may be set equal to 
unity, as in equation (39). 
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Integration of equation (39) and simplification yield 


where 


and 



(41) 


(42a) 


(42b) 


Equations (42a) and (42b) are elliptical integrals of the second kind which, when evaluated 
for a > = 90°, 120°, and 180°, yield the values of f shown in table IV. 

The values thus obtained for f^ are less than those that are obtained by including 
radiation reflected from the adjacent panel. The expression derived for f is, however, 
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TABLE IV. - CONFIGURATION FACTORS 


Dihedral angle 
between radia- 
tor panels, 

deg 

Quantity proportional 
to radiation intensity 
(eq. (39)), 

Configuration 
factor (eq. (38)), 

7r/2 

2.484 

0. 793 

2tt/3 

2.902 

.925 

TF 

3. 141 

1.0 


considered appropriate in the present case because it represents the lower bound for the 
panel arrangements presented. 
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APPENDIX C 


INERTIAL PRESSURES 


In a capillary radiator inertial pressure variations of long duration or high intensity 
can have undesirable consequences. If the radiator is designed to operate at very low 
pressures, then its incorporation into a space vehicle that is subjected to high maneuver- 
ing accelerations may be inadvisable because of possible overpressures. Excessive 
pressure at puncture sites would inevitably result in liquid escape. In general, inter- 
planetary missions are expected to entail rather low accelerations (ref. 17). Anticipated 
accelerations need to be considered in selecting a capillary dimension and header length 
to assure that pressure due to acceleration can be supported by surface tension forces. 

The greatest potential contribution to the inertial pressure is made by the longest 
liquid -filled component which is the header with length L in the radiator model de- 
scribed in the section Heat Transfer and Fluid Dynamic Considerations. The maximum 
inertial pressure Pg arising when the vehicle acceleration g occurs in the direction 
of the header axis is 

vpLs °(t) <43) 


where g Q is acceleration of gravity at the Earth's surface. Given the acceleration 
schedule of the vehicle during flight, comparison with the following equation: 


4cr(-cos g) /®o ' 


g 


pLg r 


k g 


P„ = P 
g 


(44) 


reveals whether the capillary diameter based on pump pressure through equation (1) is 
realistic. Or, combining equations (43) and (44) results in an expression for g in terms 
of the ratio P /P, that is, 

fc> 


_g__ 4q(-cos 6) f P g| 
So 5pLs o V p / 


(45) 


and stipulating a maximum tolerable value for P /P, for example 0. 1, yields a limit 

o 

for inflight accelerations. Also, as is evident from equation (44), assignment of values 
to g/g Q and Pg /P imposes restrictions upon the quantity 5gpL. Obviously, the 
selection of a low density liquid is desirable not only to reduce radiator specific weight 
but to increase the limit on acceleration. 
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APPENDIX D 


PUNCTURES IN TYPICAL CAPILLARY STRUCTURES 

Meteoroid damage to prospective capillary media should be assessed by examining 
punctures in representative liquid-filled radiator segments. The type of information 
needed in order to make this assessment is presently unavailable. 

The experiments described hereinafter indicate the damage that can be expected in 
an unfilled capillary medium of the type that may be suitable for a capillary radiator. 
These were low velocity projectile firings that yielded tentative information concerning 
the relationship between hole size, projectile diameter, and velocity. 

A typical result was obtained by firing a projectile at a corrugated steel target which 
represented a prospective capillary radiator panel segment. Two corrugated stainless - 
steel plates about 0. 038 centimeter thick were spot -welded together with their corruga- 
tions crossing at right angles. The projectile was a 0. 556 -centimeter stainless -steel 



Figure 13. - Simulated meteoroid penetration of foam metal target representing 
capillary radiator segment. 
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sphere traveling with a velocity of approximately 1524 meters per second. Examination 
of the target revealed a clean hole on the impact side with a diameter slightly larger than 
the projectile. The ragged hole on the reverse side was roughly twice the diameter of 
the projectile. 

The result of firing a 0. 556 -centimeter -diameter stainless -steel ball into a porous 
nickel target is shown in figure 13. The target was 0. 653 centimeter thick and its den- 
sity was about 20 percent of the density of solid pure nickel. Two adjacent impacts were 


' C-66-384 


(a) Projectile entrance side of target. 


[Q [ EO 

1 1 1 1 1 1 1 1 1 


C-66-385 


(bl Projectile exit side of target. 


Figure 14. - Simulated meteoroid impact of foam metal target representing 
capillary radiator segment. 





made at approximately 1524 meters per second against opposite sides of the target. The 
ball made a clean circular hole on the side facing the gun in each case. The craters 
formed on the opposite sides had a ragged, approximately conical shape, and were again 
roughly twice the projectile diameter. 

Another target tested was a monometallic sandwich consisting of three bonded layers 
of foam nickel. The two outer layers, approximately 0. 317 centimeter thick, consisted 
of 20 percent density nickel. The middle layer, approximately 0. 794 centimeter thick, 
consisted of 7 percent density nickel. The foam nickel sandwich was impacted in the 
same manner as the aforementioned targets described previously. As seen in fig- 
ure 14(a), the side facing the gun had a clean circular hole slightly larger than the 0. 556 
centimeter diameter of the projectile. The opposite side had a bowl-shaped cavity with 
a diameter of approximately twice the projectile diameter, as seen in figure 14(b). 

The hole - to projectile -diameter ratios produced by the low velocity experiments 
described in the previous paragraphs are not necessarily representative of ratios pro- 
duced by projectiles with higher velocities. In lieu of pertinent high velocity data, how- 
ever, a two -to -one ratio was adopted herein in the analysis of meteoroid puncture effects. 
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APPENDIX E 


EVALUATION OF CRITICAL PARAMETERS 

The limitations on the capillary diameter and estimates of radiator area attrition 
and liquid loss are discussed in terms of the hypothetical example and conditions de- 
scribed in the text in the section Illustrative Mission, Configuration, and Conditions. 


Limitations on Capillary Dimension 

Upper bounds on 5 are imposed by the Reynolds number, hydraulic pressure, and 
inertial forces (see the section Parametric analys is procedure) . Evaluation of equa- 
tion (19) yields a Reynolds number that is below the adopted criterion value of 100 for 

-2 

capillary diameters less than 10 meter. An upper bound is also imposed by the mini- 
mum practical pressure that can be used in the capillary system. The minimum pres- 
sure is dictated in part by the return header inlet pressure which is a function of the 
header dimensions. If a 7. 6 -centimeter -diameter return header 10 meters long is as- 
sumed, then the return header inlet and, hence, radiator inlet pressure should be of the 
order of 70 newtons per square meter since AP = 1/2 P. From equation (1) the capil- 
lary diameter corresponding to a pressure of 70 newtons per square meter is about 
10“^ meter. As discussed in appendix C, an upper bound is imposed on 6 by inertial 

forces that arise in a space vehicle. Using equation (44) and assuming that g/g ranges 

-2 -2 -S -2 ° 

from 10 to 10 yields 6„ ranging from 1. 5x10 to 1. 5x10 meter, respectively. 

& _o 

An upper bound of 6=1. 5X10 meter approximately, corresponds to the low accelera- 
tions associated with electrically propelled interplanetary vehicles (ref. 17). The pre- 

-3 

vious conditions therefore collectively dictate an upper bound on 6 between 1. 5x10 to 
-2 

10 meter. 

Given the dimensions and materials of this example, the lower bound imposed by 
conductive transfer within the panels is considerably less than that imposed by the in- 
crease of specific weight that accompanies the reduction of capillary diameter. Table I 
indicates the increase in specific weight with the decrease in diameter for titanium radi- 
ator panels. The second column indicates the number of liquid layers n^ and the first 
column the associated capillary diameter obtained by solving equation (37) for 6. The 
specific weights w. in the third column for panels of titanium are obtained with equa- 

C j 

tion (21) by assuming a surface plate thickness t_ of 2. 54x10“ meter and t /6 = 10" . 
The same value for t /6 was assumed in using equation (22) for computing the liquid 
specific weights to obtain the total specific weight of a lithium -filled panel w in the 
last column (where w = w c + w^, eq. (20)). 
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Note that the figures in table I are estimates of weight with respect to electric power 
of only the radiating panels and contained liquid and do not include the specific weights of 
the headers, manifolding, pumps, and support structures associated with the radiator. 
Note also that the specific weight of the liquid adds significantly to the total specific 
weight (compare the last two columns). The specific weights quoted in table I would be 
reduced by a substantial amount in the case of the capillary medium shown in figure 9 if 
the channel and layer partitions were discontinuous instead of solid. 

_ 4 

Upon the basis of the foregoing evaluation, 5 has a lower bound of roughly 1. 6x10 
meter if a maximum specific weight of about 5 kilograms per kilowatt (electric) is pre- 
scribed. Consequently, upon this arbitrary basis and the upper bounds mentioned pre- 
viously, 6 lies in the range from 1 . 6x10”^ to 1 . 5x10"^ meter approximately. 

Surface Area Losses and Liquid Inventory Losses 

An important factor in estimating liquid and radiator area losses is the surface 
plate thickness t . Using the method suggested in the section Puncture criteria and loss 

ir ~ 

equations , tp is equated to p thus varying with m^, the minimum mass for which a 
complete puncture occurs. 

Radiator area attrition and fluid loss for a 1000 -day exposure are shown in table III 
as a function of meteoroid mass. To calculate these quantities, the upper limit n^ in 
equations (30) and (36) is set equal to 10 grams. Masses greater than 10 grams are ig- 
nored because due to a low rate of encounter, they are unlikely to contribute significantly 
to area or fluid attrition rates. This is apparent by inspection of the last two columns 
in table in. There is no need to specify the lower limit mj because the assumed radi- 
ator surface plate thickness determines those masses producing only partial penetrations 
which are therefore ignorable according to the criterion adopted in the text. In the first 
column of table HI, the mass range for m^ from 10 to 10 - ^ gram adequately covers 
the fields of interest. To calculate values for the remaining columns, the meteoroid 
density is set to equal 1 gram per cubic centimeter (ref. 20) and the meteoroid diameters 
are calculated from equation (26). Depths of penetration into the surface plate based on 
p/d m = 20 (ref. 13) are listed in the second column. Because the surface plate thick- 
ness tp is set equal to p for each value of p appearing in the second column, all 
masses equal to or greater than each succeeding m^ produce a complete puncture. The 
third column lists the fraction of original area lost during a 1000 -day mission computed 
by equation (30). The fraction of original liquid inventory lost by vaporization computed 
by equation (36), with x = 2, is listed in the last column. Liquid loss by displacement 
is obtained by multiplying 1/4 (i. e. , l/(2x)) times the figures in the third column of 
table III. The values in table HE are rough estimates based on the assumptions discussed 
previously. 
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The selection of the radiator surface plate thickness may be based partly on the fig- 
ures appearing in table in and partly upon the estimates of depth of corrosive attack. If 
the depth of corrosive attack during the mission were less than 2. 5x10 D meter, then 

even a minimum surface plate and partition thickness t— and t , respectively, of about 

pc 

2. 5x10 meter appears practicable. A review of pertinent corrosion data indicates that 
this assumption is realistic (refs. 18, 19, and 20). However, for plate thicknesses be- 
tween about 2. 5x10“^ and 2. 5x10"^ meter, full penetration by masses ranging down to 

-9 -12 -4 

between 10 and 10 kilogram is predicted by table HI. A thickness of 2. 5X10 meter 

results in an area attrition of 0. 03 percent and combined fluid loss by evaporation and 

displacement of 0. 2 percent, approximately, according to the estimates presented in the 

tabulation. These figures support the surface plate thickness used previously to arrive 

at the specific weights quoted in table I. Note that even if the loss estimates in table III 

are too low by a factor of 10, the losses of area and fluid would still be relatively small. 

Corroborative estimates of depth of penetration in aluminum and steel made by the 

authors of references 21 and 22, respectively, indicate it is unlikely that meteoroids 

less than 10 ® gram can penetrate plates thicker than 2. 5x10"^ meter. In view of this 

and the preceding discussion, it appears that a minimum radiator surface plate thickness 

of about 2. 5x10”^ meter is adequate for most candidate materials in order to maintain a 

low puncture frequency and low liquid loss. 
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